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Abstract: The aim of our study is to investigate the possibility of habitable moons orbiting the giant 
planet HD 23079b, a Jupiter-mass planet, which follows a low-eccentricity orbit in the outer region of 
HD 23079's habitable zone. We show that HD 23079b is able to host habitable moons in prograde and 
retrograde orbits, as expected, noting that the outer stability limit for retrograde orbits is increased by 
nearly 90% compared to that of prograde orbits, a result consistent with previous generalized studies. 
For the targeted parameter space it was found that the outer stability limit for habitable moons varies 
between 0.05236 and 0.06955 AU (prograde orbits) and between 0.1023 and 0.1190 AU (retrograde 
orbits) depending on the orbital parameters of the Jupiter-type planet if a minimum mass is assumed. 
These intervals correspond to 0.306 and 0.345 (prograde orbits) and 0.583 and 0.611 (retrograde orbits) 
of the planet's Hill radius. Larger stability limits are obtained if an increased value for the planetary 

1 /3 

mass m v is considered; they are consistent with the theoretically deduced relationship of m p . Finally, 
we compare our results to the statistical formulae of Domingos et al. (2006), indicating both concurrence 
and limitations. 
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1 Introduction 



The existence of planets in orbit about solar-type stars 
is now a well-established observational result. Obvi- 
ously, the ultimate quest of this type of study is to 
discover Earth-type planets and moons located in the 
habitable zones (HZs) of their host stars. Despite a 
lack of discoveries, the existence of Earth-mass planets, 
including those hosted by solar-type stars, is strongly 
implied by various observational findings including the 
occurrence and mass distributio n of close-in super - 
Earths, Neptunes , and Jupiters ( Howard et alJl2Qlch . 
Measurements bv lHoward et al] l 2010T l indicate an in- 
creasing planet occurrence with decreasing planetary 
mass rn p akin to nip 0,48 , implying that 23% of stars 
harbor a close-in Earth- mass planet (ranging from 
0.5 to 2.0 M®); see also iMarcv fc Butlerl <|2000f ) for 
earlier results. More r ecent results were given by 
IWittenmver et all (|201lT ); they indicate based on the 
analysis of 67 solar-type stars that nearly 20% of stars 
are expected to host planets with a minimum mass of 
less than 10 M®. 

There is also a range of studies focusing on exoso- 
lar moons, especially moons possibly located in stel- 
lar HZs. Previous studies on ex os olar moo ns include 
work by IWilliams et al.1 l|l997l ). iBarnes 
(2002), and more r e cently by, e.g 
(2009), [Kalteneggerl (l2010Tl. and IWeidner fc Home 
( 2010D . IWilfiams et al.l (|1997i ) investigated appropri- 
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ate orbital parameters of possible moons, and ar- 
gued that the moons need to be large enough (i.e., 
> 0.12 M®) to retain a substantial and long-lived 
atmosphere and would also need to possess a signif- 
icant ma gnetic field to b e con tinuously habitable . Th e 
work by IKipping et all (I2009T ); iKalteneggel (|2010T ): 
Wci dner fc Horna l|2010l ) focused on the properties and 
detectability of moons of transiting planets, particu- 
larly thos e which can be detect ed as part of the Kepler 
mission. IKipping et al.1 l|2009T ) pointed out that HZ 
exomoons down to 0.2 M® may be detected and ap- 
proximately 25000 stars could be s urveyed for HZ exo - 
moons with Kepler's field of view. iKalteneggerl (|2010h 
conveyed that transmission spectroscopy of moons de- 
tected as part of the Kepler mission has a unique po- 
tential to screen t hem for habit ability in th e near fu- 
ture. In addition, Wei dner fc Horna (|2010l ) explored 
the limits on the orbits and masses of moons around 
currently known transiting exoplanets, particularly 
those discovered by the COROT mission. Recently, 
IKipping et al.l (|2013h showed seven possible candidates 
for transiting exomoons using observations from the 
Kepler mission. 

Additional recent work with implications on the 
possibility o f habit ab le exomoons w er e giv en by 
Tinnev et all J201llL iRobertson et al.1 J2012ft , and 



Heller fc Barnes! (|20ld ). iTinnev et al.1 (|201lh focused 
on the detection of a gas-giant exoplanet around 
HD 38283, a G0/G1 V star, in an eccentric orbit with 
a period of almost exactly one year; a detailed dis- 
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cussion about the potential existence of exomoons and 
exo- Trojans (see below) is conveyed by the authors. 
Another interesting study about gaseous planets about 
solar-type stars permitting the pri ncipal possibility 
of hab it able exomoons w a s give n by iRobertson et al.l 
(|2012h . iHeller fc Barnes! (|2013l ) provide a thorough 
analysis of exomoon habitability constrained by stellar 
illumination and tidal heating. 

Moon formation scenarios based on a Solar System- 
like nebula have shown that a mass scaling ratio, which 
is ~ 10 -4 , exists preventing gaseous planets using core- 
accretion to gain enough material for forming Earth- 
mass moons (Canuo & Ward 2006). However, this has 
not dissuaded the possibility of gaseous planets from 
acquiring large Earth-mass moons. Previous stud- 
ies have indicated that there can be large stability 
zones for irregular or captured satellites. In the So- 
lar S ystem, Kuiper belt objects (see, e.g., iLuu et al] 
I1997T ) have demonstrated retrograde motion result- 
ing from formation scen arios with proper conditions 
(|Schlichting fc Saril 120081 ). Moreover, many irregular 
moons have been discover ed where collisions h ave ac- 
creted to substantial sizes l|Nesvornv et al.ll2003h . Cur- 
rently, the Hunt for Exomoons with Kepler (HEK) 
(|Kipping et al.l 2012) has utilized the hypothesis of ac- 
quiring Earth-mass exomoons via capture which would 
most likely be identified as orbiting in retrograde. 

Long-term orbital stability in stellar HZs is a 
necessary condition for the evolution of any form 
of life, particularly intelligent life. Examples of 
previous studies on the orbital stability of hy- 
pothetical Earth- mass planets i n ste llar HZs in- 
clude the work bvlGehman et aD_ f 1996") 
(|2002f ). iMenou fc Tabachnild J2003ft 



..200a),lAsghari et all l|2004fl.lJones et al] 
Jones fc Sleed (l2010h. iKopparapu fc Barnesi 
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2005, 2006); 
- (2010) 

and Eber le et alJ (|201lh . These studies explore the 
dynamics of multi-body systems noting that many ex- 
trasolar planetary systems contain more than one ex- 
trasolar giant planet (EGP), which is of pivotal impor- 
tance for the orbital stability of possible Earth-mass 
object (s) in those systems. 

For example, in the HD 23079 system, a Jupiter- 
mass planet is found to orbit the star in the HZ, 
therefore jeopardizing the existence of habitable ter- 
restri al planets. Thi s system has previously been stud- 
ied bvlDvoraket al] (|2004T l. lSchwarz et alJ I^OOTT ). and 
lEberle et al] i|201ll ). The study of HD 23079 is of 
particular interest because (1) the mass and spectral 
type of the center star is fairly similar to the Sun, (2) 
the Jupiter-mass planet is found to orbit the star in 
a nearly circular orbit, and (3) any exomoon of that 
planet, if such a moon exists, would therefore exhibit 
a relatively stable climate. 

In systems like HD 23079, there is the general possi- 
bility of habitable Trojan planets as well as habitable 
exomoons. A Trojan planet is one that librates around 
one of the Lagrangian points L4 and L5 of the giant 
planet. These points lie on the giant planets orbit, 
ahead (L4) and behind (L5) the planet, each forming 
an equilatera l trian gle with the planet and its star. 
iDvorak et al] (|2004T ) investigated the stability regions 
of hypothetical terrestrial Trojans librating around L4 



and L5 in specific systems, including HD 23079, and 
obtained relationships between the size of the stability 
regions and the orbital paramet ers of the gian t plan - 
ets, particularly its eccentricity. ISchwarz et al] (| 20071 ) 
thereafter identified numerous exoplanetary systems 
that can harbor Trojan planets in stable orbits. All 
these systems also provide the opportunity for host- 
ing exomoons around Jupiter-type planets, including 
exomoons located i n stellar HZs. Upda ted results for 
HD 23079 given by lEberle et al] (|201lh showed that 
this system is also well suited for the existence of hab- 
itable Earth-type Trojans, especially if a generalized 
stellar HZ is assumed, which in the Solar System corre- 
sponds to 0.95 to 1.67 AU from the Sun feasting et al.l 
1 1993 ); see Sect. 2.2 for details. 

A more recent example of a system with the poten- 
tial of hosting both exomoons and exo- Trojans in the 
"habitability phase space" is HD 38283 as described by 
iTinnev et al] (|201ll ). This system is known to possess 
an EGP around its host star largely or entirely located 
in the stellar HZ considering that HD 38283b's orbit is 
highly eccentric (i.e., e p ~ 0.41), a notable difference to 
HD 23079b considered in the present study. Exomoons 
and exo- Trojans as companions of gas-giant exoplanets 
in 0.5-2 AU orbits (somewhat depending on the stellar 
spectral type) thus offer the possibility of providing cir- 
cumstellar habitability. Exomoons and/or exo- Trojans 
may be able to form in-situ or be captured; regard- 
ing the latter, examples incl ude the outer satellites 
of Ju piter and Saturn (e.g., iJewitt fc Haghighipourl 
120071 1. 

In the present study, we focus on the system of 
HD 23079 with a concentration on the principal possi- 
bility of habitable exomoons. Our paper is structured 
as follows. First, we discuss the stellar and planetary 
parameters of the HD 23079 system. Next, we com- 
ment on the circumstellar habitability of HD 23079. 
Thereafter, we present our results and discussion, in- 
cluding a detailed statistical analysis of the exomoon 
orbital stability limits as wel l as a comparison of ou r 
results with the formulae of iDomingos et al] |2006). 
Finally, we convey our summary and conclusions. 



2 Theoretical Approach 

2.1 Stellar and Planetary Parame- 
ters 

HD 23079 has been monitored as part of the 
Anglo-Australian P lanet Search (AAPS) program 
(JTinnev et al] [2002T I that is able to perform extraso- 
lar planet detection and measurements with a long- 
term, systematic radial velocity precision of 3 m s _1 
or better. HD 23079 was identified to host a Jupiter- 
mass planet in a relatively large and nearly circular 
orbit. HD 23079 is an inactive main-sequence star; 
iGrav et al] l|2006ft classified it as F9. 5 V (see Table 1; 
all parameters hav e their usual meaning) , an updated 
result compared to lHouk fc Cowlevl (|l975h who found 
that HD 23079 is intermediate between an F8 and GO 
star. Its stellar spectral type corresponds to a mass 
of M = 1.10 ± 0.15 M©. The stellar effective tem- 
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perature and radius are given as Tgg = 6030 ± 52 K 
and R = 1.106 ± 0.022 R , respectively (jRibas et all 
2003). Thus, HD 23079 is fairly similar to the Sun, 
though slightly hotter and slightly more massive. The 
detected planet (HD 23079b) has a minimum mass of 
m p sin i = 2.45 ±0.21 Mj. Furthermore, it has a semi- 
major axis of a p = 1.596±0.093 AU and an eccentricity 
of e p = 0.102 ± 0.031 ([Butler et al.ll2006h . correspond- 
ing to an orbital period of P = 730.6 ± 5.7 days. The 
original results bv lTinnev et al.l (|2002) indicated very 
similar planetary parameters. 

The orbital parameters of HD 23079b are relatively 
similar to those of Mars, implying that HD 23079b 
is orbiting its host star in or near the outskirts of 
the stellar HZ; see discussion below. The existence of 
HD 23079b, a planet about two and a half times more 
massive than Jupiter, makes it difficult for a terres- 
trial planet to orbit HD 23079 at a similar distance 
without being heavily affected by the g iant planet; 
see re sults from previous case s tudies bv lNoble et al.l 
(|2002T ) and lYeaeer et aD (|201ll ) who focused on the 
dynamics of HD 20782, HD 188015, and HD 210277, 
which are systems of similar dynamical settings. Con- 
cerning HD 23079, previous investigations pertaining 
to habitable t errest rial Trojan planet s were given by 
iDvorak et aD (|2004l 1 and lEberle et al.l (|201lT ). 



2.2 Stellar Habitable Zone 

The extent of the HZ of H P 23079 can be cal c ulated 
following the formali sm by | Underwood et al.l (|2003t ) 
based on the work bv lKasting et al] l|l993h . while also 
adopting an appropriate correction for the solar effec- 
tive temper ature used in the definit ion of solar lumi- 
nosity, I/©. lUnderwood et all (|2003T ) supplied a poly- 
nomial fit that allows to compute the extent of the con- 
servative and the generalized HZ based on prescribed 
values of the stellar luminosity and the stellar effec- 
tive temperature (see Table 2). Noting that HD 23079 
is considerably more luminous than the Sun, it is ex- 
pected that its HZ is more extended than the solar HZ, 
for which the inner and outer limit of the generalized 
HZ were obtained a s 0.84 and 1.67 AU, respectively 
(|Kasting et al.lll993h . 

The generalized HZ is defined as bordered by the 
runaway greenhouse effect (inner limit), where water 
vapor enhances the greenhouse effect thus leading to 
runaway surface warming, and by the maximum green- 
house effect (outer limit), where a surface temperature 
of 273 K can still be maintained by a cloud-free CO2 
atmosphere. The inner limit of the conservative HZ is 
defined by the onset of water loss, i.e., the atmosphere 
is warm enough to allow for a wet stratosphere from 
where water is gradually lost by photodissociation and 
subsequent hydrogen loss to space. Furthermore, the 
outer limit of the conservative HZ is defined by the first 
CO2 condensation attained by the onset of formation 
of CO2 clouds at a temperature of 273 K. 

For HD 23079, the limits of the conservative HZ are 
given as 1.1378 and 1.6362 AU, whereas the limits of 
generalized HZ are given as 0.989 and 1.966 AU (see 
Fig. 1). The limits of the generalized HZ are those 



employed in our numerical planetary studies[j. The 
underlying definition of habitability is based on the as- 
sumption that liquid surface water is a prerequisite for 
life, a key concept that is also the basis of ongoing and 
future searches for extrasolar habitable pla nets (e.g., 
ICatanzarite et al.ll2006l : ICockell et al.ll2009h . The nu- 
merical evaluation of these limits is based on an Earth- 
ty pe planet with a CO2/ H2O/N2 atmosphere; see work 
bv lKasting et all ll 19931 ') and subsequent studies. 

We point out that concerning the outer edge of hab- 
itability, even less conservati ve limits have been pro- 
posed in the meantime (e. g., iForget fc Pierrehumbertl 
119971 ; iMischna et"afl 12000! ) . They are based on the 
assumption of relatively thick planetary CO2 atmo- 
spheres and invoke strong backwarming that may fur- 
ther be enhanced by the presence of CO2 crystals and 
clouds. However, as these limits, which can be as large 
as 2.4 AU in case of the Sun, depend on distinct prop- 
erties of the planetary atmosphere, they are not rele- 
vant for our study. Nevertheless, we convey this type of 
limit for the sake of curiosity (see Fig. 1) , noting that it 
has properly been adjusted to 2.75 AU in consideration 
of the radiative conditions of the planetary host star, 
HD 23079. Moreover, the significance of this extreme 
limit has meanwhile been chall enged based on det ailed 
radiative transfer simulations l|Halevv et al.ll2009l ). 

A more thorough investigation of exoplanet and ex- 
omoon habitability points to the consideration of ad- 
ditional factors that might influence the facilitation of 
an habitable environ ment . Detailed reviews by, e.g., 
lLammer et"afl (|2009f ) and iHorner fc Jonesl (|2010T ) in- 
dicate the importance of magnetic fields and plate tec- 
tonics as well as the prospects for the object of be- 
ing hydrated as potentially decisive. Exomoons of 
HD 23079b are expected to be in a fortunate situation 
from that point of view as one would expect that tidal 
effects of the EGP on the moon might result in a hot in- 
terior, and hence make plate tectonics and a magnetic 
field more likely. Moreover, if the giant plan et had mi- 
grate d in from beyond the snow line (e.g., lArmitagg 
I2007T ). it could help hydrate the moon. 

As part of our evaluation of habitability around 
HD 23079, we explored whether or not the orbit 
of HD 23079b is indeed fully located in the HZ of 
HD 23079. Thus, we calculated the limits of the gen- 
eralized HZ by also considering the uncertainties in 
the stellar luminosity and ef fective temperatur e while 
employing the formalism by lUnderwood et aD (J2003) 
based on second-order polynomials. The stellar lu- 
minosity and effective temperature were assumed as 
1.45±0.08 L Q and 6030±52 K, respectively. Therefore, 
the outer limit of the generalized HZ, denoted as HZ-o 
(gen.), see Table 2, is obtained as 1.966 ± 0.054 AU. 

In addition, we calculated the periastron and apas- 
tron of HD 23079b, denoted as a pcr and a ap , respec- 
tively, while also considering the statistical uncertain- 
ties in the planetary orbital parameters, i.e., the semi- 
major axis and eccentricity, given as 1.596 ± 0.093 AU 
and 0.102 ± 0.031, respectively (see Table 1). Hence, 



x The physical limits of habitability are much less strin- 
gent than implied by the numerical precision of these val- 
ues; nevertheless, these values were used for checking if the 
Earth-mass planet has left the stellar HZ. 
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the values for the periastron and apastron are identi- 
fied as 1.432±0.094 AU and 1.758±0.113, respectively 
(see Table 2). Our analysis shows that the orbit of 
the EGP is almost certainly completely positioned in 
HD 23079's HZ; the probability for HD 23079b's orbit 
to be partially outside the generalized HZ is only 5%. 
It should be also n oted that following the study by 
IWilliams fc Pollard! l|2002l ) brief excursions from the 
stellar HZ do not necessarily nullify planetary habit- 
ability because the latter is expected to mainly depend 
on the average stellar flux received over an entire orbit. 



3 Results and Discussion 

3.1 Orbital Stability Simulations 

For our simulations of the HD 23079 system, we con- 
sider both the observed giant planet HD 23079b and 
an object of one Earth-mass, i.e., 3.005 x 10 -6 M@, 
which based on our choice of initial conditions is ex- 
pected to constitute a moon in regard to the giant 
planet. Consequently, we execute a grid of model sim- 
ulations by considering different starting conditions for 
the moon candidate. The method of integration as- 
sumes a fourt h-order Runge-Kutta integration scheme 
(JGarcial 120001 ) and an eighth-o rder Gragg-Bulirsch - 
Stoer (GBS) integration scheme (jGrazier et al.l ll996). 
The code has been extensively tested against known 
analytical solutions, including the two-body and re- 
stricted three-body problem (see |Noble et all 1 20021 : 
ICuntz et ail 120071 ; lEberle et all 120081 . for detailed re- 
sults) . 

We limit our computations to a simulation time of 
10 6 yrs and apply a time-step of 10 -4 yrs for the in- 
tegration scheme; the latter is found to be fully ap- 
propriate. In conjunction with previous investigations 
of habitable Trojan planets in the HD 23079 system 
(|Eberle et al.ll2011l ). we pursued test simulations com- 
paring the planetary orbits based on three different 
integration time-steps, which were: 10 -3 , 10 -4 and 
10~ 5 yrs. In particular, we evaluated ARij, i.e., the 
magnitude of the difference between the position of 
the planet when different step sizes of 10 _I and 10 - -' 
were used. We found no discernible differences in the 
outcome between models with time-steps of 10 -4 and 
I0 -5 yrs, thus justifying the usage of 10 -4 as time-step 
for this type of simulations. 

The giant planet always starts out at the 3 o'clock 
position relative to the star placed at the origin, which 
is assumed to coincide with its periastron position. 
The initial condition (i.e., starting velocity) for the 
orbit of the terrestrial moon is chosen such that it is 
assumed to be coplanar regarding the motions of the 
star and the EGP, as well as in a circular orbit about 
the EGP, although it is evident that it will be signif- 
icantly affected immediately by gravitational pull of 
the star, which is expected to prevent from continu- 
ing its circular motion. Furthermore, the moon always 
starts at the 9 o'clock position relative to the giant 
planet. Our set of models is defined by five values for 
the semimajor axis a p and five values for the eccen- 
tricity e v of the giant planet, which range from 1.503 



to 1.689 AU in increments of 0.0465 and from 0.071 
to 0.133 in increments of 0.0155, respectively. For the 
EGP, we assume its minimum mass value of 2.45 Mj 
as default. However, for models with a p = 1.596 AU 
and e p — 0.102, the mass of the EGP is increased by 
factors of 1.5 and 2.0. 



3.2 Derivation of Orbital Stability 
Limits 

In our model simulations we consider both prograde 
and retrograde orbits of the moon candidate. In the 
selected coordinate system the giant planet orbits the 
star in a counter-clockwise sense, which means that for 
prograde orbits, the moon orbits the giant planet in a 
counter-clockwise sense as well. Hence, for retrograde 
orbits, the moon orbits the giant planet in a clockwise 
sense. The overwhelming majority of our simulations 
assumes a minimum mass for the giant planet; these 
models will be the focus of the following methodology. 

The identification of model-dependent orbital stabil- 
ity limits for the moon requires a large number of simu- 
lations since the stability limits (i.e., maximum radii of 
possible orbits) can only be determined through trial- 
and-error. As the selected approach, both for prograde 
and retrograde orbits, a two-step process is adopted. 
First, we consider a preliminary 3x3 grid for the tar- 
geted array of values for the semimajor axis a p and 
eccentricity e p of the EGP; for a p , they are given as 
1.503, 1.596, and 1.689 AU, whereas for e p , they are 
given as 0.071, 0.102, and 0.133. 

For each grid value, a simulation time of 5 x 10 4 yrs 
is used to obtain estimates for the stability limit of 
each (a p , e p ) combination. For each prograde model, 
we perform a total of 201 simulations, based on in- 
tervals for the moon's starting distance with 0.0520 
and 0.0720 AU as lower and upper limits, respectively. 
Thus, the simulations are advanced in increments of 
1 x 10~ 4 AU regarding the moons starting distance. 
For each retrograde model, we simulate a total of 201 
initial conditions as well; here the limits of the intervals 
for the moon's starting distance are given as 0.1050 and 
0.1250 AU. This approach amounts to a total of nearly 
3600 simulations for the prograde and retrograde mod- 
els combined. 

To refine our search for orbital stability limits, a sec- 
ond step is implemented involving two aspects, which 
are: augmenting the preliminary coarse 3x3 grid to 
the full 5x5 grid, while also considering substantially 
longer simulation times for each model, i.e., target 
timespans of 10 6 years. For each (a p , e p ) combination, 
these long-term simulations again utilize intervals for 
setting the moon's starting distance. For each grid 
combination, we execute a total of 16 simulations cen- 
tered at the expected value for the model-dependent 
orbital stability limit. 

For models already part of the preliminary 3x3 
grid, the expected values for the orbital stability lim- 
its are obviously given by the results of the previous 
short-term computations. For models beyond the pre- 
liminary 3x3 grid, the expected values for models with 
intermediate values of semimajor axis (i.e., 1.5495 and 
1.6425 AU) are generated through applying a poly- 
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nomial fitting assuming constant values for a p (i.e., 
1.503, 1.596, and 1.689 AU). The polynomial fitting 
process also assisted us in estimating where stability 
limits should occur for intermediate values of eccen- 
tricity (i.e., 0.0865 and 0.1175). 

Our results are given in Tables 3, 4, and 5; for con- 
venience, the data of Table 3 and 4 are also depicted 
in Figs. 2 and 3. Based on the total number of simula- 
tions, it is found that, if the EGP is assumed to possess 
minimum mass, the stability limit for habitable moons 
varies between 0.05236 and 0.06955 AU (prograde or- 
bits) and 0.1023 and 0.1190 AU (retrograde orbits) de- 
pending on the orbital parameters of the giant planet. 
By keeping the EGP's mass and the semimajor axis 
fixed, expr essions are encountered str ucturally similar 
to those of iHolman fe Wiegertl (11999T ); however, with 
different coefficients naturally related to our particular 
study. 

Somewhat larger stability limits of the candidate 
moon are obtained if an increased value for the mass 
of the EGP is considered. Specifically, if the mass of 
the giant planet is increased by factors of 1.5 and 2.0, 
the stability limit for prograde orbits increases from 
0.05712 to 0.06290 and 0.06940 AU, respectively. Fur- 
thermore, for retrograde orbits, the stability limit of 
habitable moons is found to increase from 0.1109 to 
0.1291 and 0.1400 AU, respectively. 

3.3 Examples 

Examples of orbital stability and instability of exo- 
moons around HD 23079b are given in Figs. 4 to 
8. The orbital parameters of the giant planet con- 
sidered in these case studies are a p = 1.596 AU and 
e p — 0.102; additionally, the EGP is assumed to be of 
minimum mass. Figures 4 and 5 indicate that for pro- 
grade orbits short-term (~ 10 4 years) stability — sub- 
sequently confirmed to also constitute long-term sta- 
bility; see Sect. 3.2 — is found for a separation distance 
of 0.0574 AU, and instability is found for a separation 
distance of 0.0576 AU. The escape occurs after only 
4442 years. Furthermore, as indicated in Figs. 6 and 
7, short-term stability is found for retrograde orbits for 
a separation distance of 0.1112 AU, whereas instability 
is found for a separation distance of 0.1113 AU. Here 
the escape occurs after 3145 years. Figure 8 shows 
an example of a long-term (~ 10 6 years) instability 
transition for an exomoon in a retrograde orbit. Here 
we chose two nearby starting positions, i.e., 0.1107 AU 
and 0.1111 AU, to determine if the choice of initial con- 
ditions in this region will largely determine the long- 
term stability. 

Figure 8 (top panel) depicts the orbit, within a rotat- 
ing reference frame, of an exomoon in a retrograde or- 
bit using the 0.1111 AU starting position; the exomoon 
escapes from the system after 4000 years. Further- 
more, we investigated how the semimajor axis (mid- 
dle panel) and eccentricity (bottom panel) evolved over 
time considering these two nearby initial positions. 
Figure 8 (middle panel) shows that the semimajor axis 
for a m — 0.1107 AU (depicted as dots) remains rela- 
tively constant with minor fluctuations in time for at 
least 10 5 years, whereas the alternative a m (depicted 



as plus symbols) experiences an increase. The subse- 
quent decrease in semimajor axis indicates the loss of 
the moon due to ejection as also seen in the evolution 
of eccentricity. 

Figure 8 (bottom panel) demonstrates more clearly 
for a 10 4 year timespan the instability transition oc- 
curring at 4000 years for the a m = 0.1111 AU case, 
exhibiting a similar trend in the value of eccentric- 
ity. These simulations were performed with the thor- 
oughly tested GBS integration scheme. Incidentally, 
in o ur previous study a imed at habitable Trojan plan- 
ets l|Eberle et al.ll201lh . we already encountered a by- 
chance scenario with a habitable moon in a retrograde 
orbit about the giant planet; it occurred due to a par- 
ticular set-up of the initial conditions. In this case, 
the semimajor axis of the exomoon was identified as 
0.051 AU, which is consistent with the stability limit 
attained in our present, more detailed study. 

Technically speaking, the results on limits of orbital 
stability described so far should be referred to as outer 
orbital stability limits. The inner orbital stability lim- 
its of exomoons in the HD 23079 system are given by 
the Roch e limit of the giant planet; see discussion by 
IWilliama (|2003l ) and references therein. The evalua- 
tion of the Roche limit typically yields a distance of 
about 2.45 radii from the central object as stability 
limit; its assessment requires an approach distinctly 
separate from the main focus of this paper. The Roche 
limit is known to depend on, e.g., the radius, oblate- 
ness, and density of the primary, and structural and 
rotational properties of the satellite. In the absence 
of a more detailed discussion, the inner stability limit 
for exomoons around HD 23079b can be estimated as 
about 0.0015 AU; note that this estimate is weakly 
dependent on the mass of the giant planet. 

3.4 Statistical Assessments 

Another aspect of our investigation concerns the statis- 
tical assessment of the orbital stability limits obtained 
through numerical simulations. Previously, this has 
been done by applying the c r iterion of Hill s tability 
(e.g- lHamilton fe Burnslll992l ; lDo"nnisonll2010h , which 
often serves as a customary approximation utilized in 
conjunction with time-dependent simulations to gauge 
the onset of or bital instability of planets i n multi-body 
systems fe.g.. [Menou fe Tabachnikl 120031 : iJones et al.l 
120051 : Ivon Bloh et al.1 120071 ). although the Hill stabil- 
ity concept has previously been ide ntified as largely 
defunct in particular applications l|Cuntz fe Yeagen 
120091 ). 

For the parameter space as studied (see Tables 3 
and 4) we find that the stability limit for the exo- 
moons varies between 0.05236 and 0.06955 AU (pro- 
grade orbits) and between 0.1023 and 0.1190 AU (ret- 
rograde orbits) depending on the orbital parameters of 
the Jupiter- type planet if a minimum mass is assumed. 
In general, we find that the larger the semimajor axis 
a p , the larger the stability limit. Moreover, we also 
find that the larger the eccentricity e p , the smaller the 
stability limit. 

Guided by the definition of the Hill radius (see Sect. 
3.5), we assume that the exomoon's orbital stability 
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limit ziim follows a relationship akin to 

Zlta = ^(g^l-e^ , (1) 

with a p , e p and m p as defined. The parameters Ao, a, 
13, and 7 are to be determined from the data obtained 
by our orbital stability simulations. 

First we focus on the case where the EGP is as- 
sumed to have minimum mass, m p ; this allows us to 
constrain all parameters but 7 (see Table 6). In this 
case, there are 25 models of prograde orbits and 25 
models of retrograde orbits available. Fitting the data 
yields A = 0.296, a = 1.02, and /3 = 0.68 for the pro- 
grade models, and A = 0.572, a = 0.99, and /3 = 0.52 
for the retrograde models. As quality checks we as- 
sessed the mean and maximum errors obtained for the 
two data sets as well as the errors for the master mod- 
els (MMod) given as a p = 1.596 AU and e p = 0.102. 

We find that these fits for the prograde and retro- 
grade cases are fully acceptable; the errors are typically 
< 10% while noting that in the retrograde case the er- 
rors are considerably smaller as some stability limits in 
the prograde case are adversely affected by resonance 
phenomena (see below). In order to simplify the fits, 
we pursued tests with a = 1.0 and ft — 0.5 as well as 
f3 — 1.0 for the prograde and retrograde cases. Even 
though these choices were found to be largely inconse- 
quential for the quality of the fits, we elected a — 1.0 
and /3 = 1.0 (adjusted fit 2) as our preferred choice 
for both the prograde and retrograde cases, as it also 
reflects the functional form of the Hill stability radius 
(see Sect. 3.5). Additionally, it is found that the pre- 
ferred fit renders the lowest maximum and MMod er- 
rors compared to the best fit and the adjusted fit 1. 
Moreover, the differences regarding the mean errors 
between adjusted fit 2 (our preferred choice) and the 
best fit are also found to be relatively minor. 

Some values in Table 3 do not perfectly fit the 
expected progression of stability limit regarding the 
surrounding (a p , e p ) results. Specifically, we find 
that for the cases where a p — 1.6425, the predicted 
stability limits would reside near an 11:2 mean mo- 
tion resonance. Recent studies with the circumbi- 
nary planet, Kepler-16b, has shown that the Satur- 
nian planet routinely clears the s urrounding stability 
pockets near the same resona nce dQuarles et al.ll2012l ; 
iPopova fc Shevchenko | [ 20 12!). Thus our estimations 
of the stability limit for this region would be trun- 
cated due to this type of behavior. Fortunately, our 
overall statistics is not considerably impacted by this 
phenomenon. 

As part of our study, the mass of the EGP was in- 
creased by factors of 1.5 and 2.0 relative to its mini- 
mum mass (see Table 5). This also allows us to de- 
duce the parameter 7 (see Eq. 1), which we identified 
as 7 = 0.309 ± 0.061 (la) with a M = 0.031, the stan- 
dard error of the mean. Note that the value of 7 as 
determined does not significantly deviate from the the- 
oretical value of 7 = 1/3; it obviously agrees with the 
model-dependent value by less than Ictm- 

Another quantity of interest concerning our stability 
studies is the ratio of orbital stability limits between 
models of moons in retrograde and prograde orbits, 



retrograde / prograde 



n = «lTm"° ~~"IZnm'" ■ ( 2 ) 

We find that n is given as 1.891±0.094 for models with 
EGP's minimum mass; here ctm is given as 0.019. If 
the four models with the increased mass of the EGP 
are also included, n is identified as 1.902 ± 0.098 with 
ctm = 0.019. Thus, in conclusion, the stability limit 
for exomoons in retrograde orbits is increased by 90% 
over the stability limit for prograde orbits. 

3.5 Comparison to Previous Work 

Finally, we can also compare the stability limits ob- 
tained through our numerical simul ations to the re- 
sults of the fitting formulae given bv lDomingos et aU 
(2006) (D06). With R H denoting the wel l-established 
Hill radius (e.g.. Hamilton fc BurnsHl992h . given as 



Bu 



<m 



1/3 



(l-e p ) 



(3) 



where a p and e p are the semi-major axis and eccentric- 
ity of the planetary orbit, respectively, and M and m p 
are the masses of the star and planet, respectively, the 
limits of stability for prograde and retrograde orbits of 
the Earth-mass exomoon are given as 

prograde = q.4895 (1.0-1.0305%, -0.2738e m )-R H (4) 

and 



retrograde 
2 lim 



0.9309 (1.0 - 1.0764e p - 0.9812e„ 
+0.9446e p e m )i? H , 



(5) 



respectively, with e m denoting the orbital eccentricity 
of the exomoon. 

A simple test of using the D06 formulae for the 
limiting case of e p — > and e m — > reveals a 
promising result: the value of n (D06) is identified 
as 0.9309/0.4895 = 1.902, which is in perfect agree- 
ment with the n value identified for our simulations 
(see Sect. 3.4). 

However, a more appropriate evaluation of the D06 
formulation requires the checking of its applicability to 
the Ziim values for prograde and retrograde orbits in a 
separate manner. Thus, we introduce 



D06 / this work 



/ z llm 



(6) 



In this regard we find a significantly different outcome 
for prograde and retrograde models. If e m = is as- 
sumed, no agreements between our numerical simula- 
tions and the results of the D06 formulae are attained. 
However, during orbital instability transitions, the or- 
bital eccentricity of the exomoon is expected to be rel- 
atively large, i.e . . e m — > 1; see, e.g., the analysis by 
lEberle fc Cuntj (|2010al ). 

Therefore, for prograde orbits, if e m = 0.95 is used, 
we find q = 1.027 with ct m = 9.11 x 10" 3 . More- 
over, if e m = 0.99 is used, we find q — 1.005 with 
ctm = 9.06 x 10~ 3 . Hence, for prograde orbits, the 
D06 formula is identified as fully applicable to repre- 
sent the Zy lm values of orbital stability for exomoons 
in the HD 23079 system. However, similar tests for 
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retrograde orbits show that the DOB formula is off by 
about an order of magnitude. Thus, we conclude that 
the D06 formulation has not been well-tested for pa- 
rameter sets of the latter kind. 



4 Summary and Conclusions 

The aim of our study was to explore the general possi- 
bility of habitable moons in the HD 23079 star-planet 
system. This system consists of a main-sequence star 
slightly hotter than the Sun. It also contains a Jupiter- 
mass planet, HD 23079b, with a minimum mass of 
2.45 Mj in a slightly elliptical orbit that is almost cer- 
tainly completely positioned in HD 23079's HZ. We 
studied if Earth-mass habitable moons can potentially 
exist in this system. Thus, we deduced the orbital 
stability limits of moons orbiting HD 23079b by con- 
sidering both prograde and retrograde orbits. 

The set of models of our study also took into ac- 
count the observational uncertainties of a p and e p for 
HD 23079b, which for the minimum mass case of the 
EGP resulted in the establishment of a 5 x 5 grid of 
stability limits for both prograde and retrograde or- 
bits; clearly, a much larger number of models was pur- 
sued for the derivation of the stability limits due to 
the required test simulations given by different exo- 
moon starting distances from the EGP. Moreover, the 
semimajor axis a p of the EGP was varied from 1.503 
to 1.689 AU in 0.0465 increments, whereas its eccen- 
tricity e p was varied from 0.071 to 0.133 in 0.0155 
increments. Furthermore, for the master model with 
a p = 1.596 AU and e p = 0.102, the mass of the EGP, 
m p , was increased by factors of 1.5 and 2.0. 

For the minimum mass case of the EGP, the orbital 
stability limits were found to vary between 0.05236 
and 0.06955 AU (prograde orbits) and 0.1023 and 
0.1190 AU (retrograde orbits), corresponding 0.306 
and 0.345 (prograde orbits) and 0.583 and 0.611 (ret- 
rograde orbits) of the EGP's Hill radius. Somewhat 
larger stability limits for the exomoons were obtained 
if an increased value for the mass of the Jupiter-type 
planet was adopted, as expected. Generally, it was 
found that the larger the semimajor axis a p , the larger 
the stability limit. It was also found that the larger the 
eccentricity e p , the smaller the stability limit. A small 
number of exceptions was identified when the orbital 
motion of the exomoon was affected by resonances. 
Also note that the limits as deduced constitute, tech- 
nically speaking, outer orbital stability limits; the in- 
ner limits of orbital stability are, in essence, given by 
the EGP Roche limit. 

Another element of our study was the derivation 
of the functional dependencies of the orbital stabil- 
ity limits on the EGP parameters based on zy lrn oc 
a p m p (l — e p y guided by the customary notion of Hill 
stability. Based on our simulations, we adopted a — 1, 
/3 = 1, and 7 = 1/3 for both prograde and retrograde 
orbits, which is in agreement with the functional de- 
pendencies of the Hill stability radius. In addition, 
we also found that the stability limits for prograde 
orbits wer e in perfect agre e ment with the analytical 
formula of iDomingos et all ([2006). However, no such 



agreement was identified for retrograde orbits, which 
is most likely attributable to that the formula fo r ret- 
rograde satellite motion of lDomingos et ahl (J2006) may 
be inapplicable for e p — ► and e m — > 1. 

A further quantity considered in our study concerns 
the ratio of orbital stability limits between models of 
moons in retrograde and prograde orbits r\ (see Eq. 2). 
According to our model simulations, we found that 77 
is given as 1.902 with ayi = 0.019. Thus, on average, 
the stability limit for retrograde orbits is increased by 
90% compared to the limit of stability for prograde or- 
bits. This result is, in pr i nciple , consistent with the 
previous study bv ljefferva IJ1974T ) who explored the or- 
bital stability of the restricted 3-body problem based 
on the usage of Henon stability maps. He identified 
considerably broader stability regions for the massless 
particle-type object if a retrograde instead of a pro- 
grade orbit about the primary object was assumed. 

Another example for significantly increased stabil- 
ity limits for r e trograd e orbits has been given by 
lEberle fc Cunt d ((2010b). They pursued a theoretical 
investigation of a previous observ ational finding about 
the v Octantis binary system IJRamm et al.l [2009) 
that indicates the possible existence of a Jupiter-mass 
planet, although the planet appears to be located out- 
side the zone of orbital stability. Based o n detailed nu- 
merical studies, lEberle fc Cunta (|2010bh argued that 
the Jupiter-mass may be in a retrograde orbit relative 
to the motion of the binary components, which would 
allow the planet to be orbitally stable. Additional re- 
sults in line with t his reasoning hav e recently been 
given by Quar les et alj (|2012| ) and iGozdziewski et aU 
(|2013h . 

The study of exomoons is currently receiving height- 
ened attention because of the ongoin g Kepler mis- 
sion. For example, recent work by iKipping et alj 
(2009) addressed the detectability of moons of tran- 
siting planets, including those that are expected to be 
detected by Kepler or photometry of approximately 
equal quality. These studies convey both the pre- 
dicted transit timing signal amplitudes and the es- 
timated uncertainties on such measurements to pro- 
vide realistic confidence levels for the detection of 
such bodies for a broad spectrum of orbital arrange- 
ments. They concluded that HZ exomoons down to 
0.2 Mq, may be detectable and approximately 25000 
stars could be surveyed for H Z exomoons wi t h Ke - 
pler's field of view. Moreover, IKipping et al.1 (1 20131 ) 
recently presented seven possible candidates for tran- 
siting exomoons using observations from the Kepler 
mission. Previously, iKalteneggerl (2010) discussed the 
possibility of screening the atmospheres of exomoons 
for habitability based on the concept of biomark- 
ers. These investigations will potentially allow iden- 
tifying the existence of habitable exomoons around 
different types of stars, which will allow contest- 
ing results an d pred i ctions from earlier studi es by 
IWilliams et"al] (1199711. iBarnes fc O'Brienl (I2002T). and 
more re cently by Kipping et al. (200 9J'). | Kaltenegger 
201C h . IWeidner fc Hornd (I2010T ). iRobertson et al. 



2013 ) , iHeller fc Barnes! (|2013h ~and others 



This study complements a previous investigat ion by 
our research group, given by lEberle et al.1 (|201ll ). that 
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Figure 1: Extent of the HZ for HD 23079, denned 
by its conservative limits (dark gray) and general- 
ized limits (medium gray). In addition, we de- 
pict the outer limit of an extreme version of the 
generalized HZ (ligh t gray) following the work by 
Mischna et al.1 ( 2000f ) , although this limit may still 
be controversial. The orbit of HD 23079b is de- 
picted by a thick solid line. 



showed that the HD 23079 star-planet system is also a 
suitable candidate for hosting habitable Trojan plan- 
ets. 
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Figure 5: Separation distance of the moon from 
HD 23079b for the orbital stability study depicted 
in Fig. 4. Note that the moon has been placed in 
a prograde orbit about the planet with its starting 
distance given as 0.0574 AU (top) and 0.0576 AU 
(bottom). 
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Figure 2: Depiction of orbital stability limits of habitable moons for prograde orbits; see Table 3 for data 
information. The data are grouped in regard to the EGP's eccentricity e p . 
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Figure 3: Depiction of orbital stability limits of habitable moons for retrograde orbits; see Table 4 for 
data information. The data are grouped in regard to the EGP's eccentricity e p . 
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Figure 4: Orbital stability simulations for the moon assumed to orbit HD 23079b (red line) for an 
elapsed simulation time of 10 years. The orbital parameters of the planet are given as a p = 1.596 AU 
and e p = 0.102. The moon has been placed in a prograde orbit about the planet with its starting distance 
given as 0.0574 AU (top) and 0.0576 AU (bottom). 
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Table 1: Stellar and Planetary Parameters 



Parameter 



Value 



Reference 



Spectral Type 

RA 

DEC 

Teff (K) 

R(R & ) 

M (M ) 
M v 

A/boi 

Distance (pc) 
rripsini (Mj) 
P (days) 
a P (AU) 



F9.5 V 

3 h 39 m 43.0952 s 

- 52° 54' 57.017" 

6030 ± 52 

1.106 ± 0.022 

1.10 ± 0.15 

4.42 ± 0.05 

4.25 ± 0.05 

34.60 ± 0.67 

2.45 ± 0.21 

730.6 ± 5.7 

1.596 ± 0.093 

0.102 ± 0.031 



Gray et al. (2006) 

ESA (1997) a > b 

ESA (199JJ a ' b 

Rib as et al. (2003) 

Ribas et al. (2003 



ESA (1997)°' M 
ESA (1997)°' M 
ESA (1997) a ' M 
Butler et al. C2006) 
Butler et al. (]>006) 
Butler et al. C20061 
Butler et al. C2006 



"data from SIMBAD, see http://simbad.u-strasbg.fr 

b adopted from the Hipparcos catalogue 

c based on spectral type 

d based on parallax 28.90 ± 0.56 mas 



Table 2: Orbit of the EGP and HZ 



Parameter 


Value 


a 


flpcr 


1.432 


0.094 


flap 


1.758 


0.113 


HZ-i (gen.) 


0.989 


0.027 


HZ-i (cons.) 


1.138 


0.031 


HZ-o (cons.) 


1.636 


0.044 


HZ-o (gen.) 


1.966 


0.054 



Note: All data are given in units of AU. 



Table 3: Stability Limits of Habitable Moons; Prograde Orbits 



1.5030 1.5495 1.5960 1.6425 1.6890 

(AU) (AU) (AU) (AU) (AU) 



0.0710 


0.06186 


0.06384 


0.06572 


0.06750 


0.06955 


0.0865 


0.05410 


0.05590 


0.05760 


0.05930 


0.06110 


0.1020 


0.05376 


0.05529 


0.05712 


0.05870 


0.06035 


0.1175 


0.05241 


0.05451 


0.05591 


0.05820 


0.05900 


0.1330 


0.05236 


0.05399 


0.05562 


0.05760 


0.05815 



Note: a p is given in units of AU. The stability limits represent maximum radii for the orbits. 
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Tabic 4: Stability Limits of Habitable Moons; Retrograde Orbits 



^ p 






'h 








1.5030 


1.5495 


1.5960 


1.6425 


1.6890 




(AU) 


(AU) 


(AU) 


(AU) 


(AU) 


0.0710 


0.1060 


0.1092 


0.1125 


0.1157 


0.1190 


0.0865 


0.1051 


0.1081 


0.1114 


0.1147 


0.1179 


0.1020 


0.1041 


0.1073 


0.1109 


0.1137 


0.1169 


0.1175 


0.1030 


0.1062 


0.1092 


0.1119 


0.1156 


0.1330 


0.1023 


0.1056 


0.1086 


0.1067 


0.1099 



Note: a p is given in units of AU. The stability limits represent maximum radii for the orbits. 



Table 5: Stability Limits of Habitable Moons for Different EGP Masses 



Mass of EGP Prograde Orbit Retrograde Orbit 



1.0 x m p 
1.5 x m p 
2.0 x m 



V 



0.05712 
0.06290 
0.06940 



0.1109 
0.1291 
0.1400 



Note: The simulations are given for a p = 1.596 AU and e p = 0.102. All data are given in units of AU. 
The stability limits represent maximum radii for the orbits. 



Table 6: Functional Parameters of zu m , Error Analysis 



Model 


Type 


A 


a 


13 


Mean Error 

% 


Max. Error 

% 


MMod Error 

% 


Prograde Orbits 


best fit 


0.296 


1.02 


0.68 


2.5 


11.4 


3.2 


Prograde Orbits 


adjusted fit 1 


0.293 


1.00 


0.50 


2.8 


11.9 


3.7 


Prograde Orbits 


adjusted fit 2 


0.309 


1.00 


1.00 


2.5 


10.4 


1.9 


Retrograde Orbits 


best fit 


0.572 


0.99 


0.52 


0.5 


4.6 


3.5 


Retrograde Orbits 


adjusted fit 1 


0.568 


1.00 


0.50 


0.5 


4.7 


3.6 


Retrograde Orbits 


adjusted fit 2 


0.600 


1.00 


1.00 


1.1 


3.0 


1.9 
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Figure 6: Orbital stability simulations for the 
moon assumed to orbit HD 23079b (red line) for 
an elapsed simulation time of 10 4 years. The or- 
bital parameters of the planet are given as a p — 
1.596 AU and e p — 0.102. The moon has been 
placed in a retrograde orbit about the planet with 
its starting distance given as 0.1112 AU (top) and 
0.1113 AU (bottom). 
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Figure 7: Separation distance of the moon from 
HD 23079b for the orbital stability studies de- 
picted in Fig. 6. Note that the moon has been 
placed in a retrograde orbit about the planet with 
its starting distance given as 0.1112 AU (top) and 
0.1113 AU (bottom). 
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Figure 8: Orbital stability simulations for the 
moon assumed to orbit HD 23079b (small five-side 
star) . The orbital parameters of the EGP are given 
as dp = 1.596 AU and e p = 0.102. The moon has 
been placed in a retrograde orbit within a rotating 
reference frame about the planet with its starting 
distance given as 0.1111 AU (top). Time series 
of the semimajor axis (middle) and eccentricity 
(bottom) are shown with respect to the designated 
starting distance in a retrograde orbit. 
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